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As  part  of  a continuing  study  of  material  models  for  ground  shock  calculations 
some  aspects  of  both  rate-dependent  and  rate-independent  models  for  tensile 
behavior  are  examined.  In  particular,  the  SRI  brittle  fracture  model  is 
studied  and  compared  with  a simple  tension  cut  off  procedure.  For  large-scale 
problems/L(^octibod  herain)^  it  is  important  to  avoid  rate-dependent  brittle 
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20.  ABSTRACT  (Continued) 

the  other  hand,  rate-independent  brittle  fracture  models  are  incompatible 
with  current  continuum  uniqueness-stability  theory  and  hence  may  not  be  use- 
able  for  these  problems.  It  is  necessary  to  examine  this  situation  if  tensile 
modeling  techniques  are  to  improve.  Some  ideas  for  resolving  this  conflict 
are  discussed  and  recommendations  for  further  study  are  made. 
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INTRODUCTION 


In  seeking  to  Improve  the  tensile  modeling  of  geological  materials 
for  use  in  ground  shock  applications,  a review.  Ref.  [1],  was  made  of 
current  practices  in  the  field.  Both  rate-dependent  and  rate- 
independent  models  were  summarized  and  it  was  recommended  that  some  of 
the  more  promising  models  be  Investigated  in  somewhat  greater  detail. 

In  this  report,  the  Stanford  Research  Institute  (SRI)  brittle 
fracture  model  (rate-dependent)  is  studied  and,  to  some  extent,  is 
compared  with  a simple  tension  cut-off  procedure  (rate-independent). 
Although  the  SRI  model  was  originally  designed  for  metals,  it  has  been 
adapted  for  ground  shock  applications  in  rocks  in  recent  years.  A copy 
of  the  model  routine  BFRACT  was  obtained  from  SRI.  It  has  been 
examined  theoretically  and  computationally  for  one-dimensional  problems 
and  some  of  the  results  are  presented  herein.  The  model  was  first 
exercised  to  simulate  simple  uniaxial  extension  and  then  Incorporated 
into  the  one-dimensional  wave  code  WONDY  for  plate  slap  calculations. 

As  written,  the  BFRACT  algorithm  requires  the  storage  of  more  than 
ten  memory  parameters  per  mesh  point.  This  added  computer  storage 
requirement  is  a major  impediment  to  its  direct  incorporation  into 
large-scale  ground  shock  calculatons.  In  Section  II,  modifications  to 
the  algorithm  are  suggested  to  ameliorate  this  situation,  following  a 
review  of  the  continuum  assumptions  upon  which  BFRACT  is  based.  In 
Section  III,  one-dimensional  numerical  tests  of  BFRACT,  wlch  and  without 
modifications,  are  presented.  It  is  shown  that  BFRACT  is  highly  rate- 
dependent,  but  much  less  sensitive  to  the  variation  of  many  of  tlie 
BFRACT  material  parameters. 


In  Section  IV,  it  is  pointed  out  that  time  steps  greater  than  100 


microseconds  arc  used  In  typical  large-scale  ground  sliock  calculations, 
but  that  time  scales  of  less  than  10  microseconds  are  believed  necessary 
to  resolve  the  rate-dependent  tensile  behavior  of  geological  materials. 
For  practical  computations,  therefore,  it  is  important  to  avoid  rate- 
dependent  models,  if  possible.  On  the  other  hand,  it  is  pointed  out  that 
current  continuum  uniqueness-stability  theory  is  incompatible  with  rate- 
independent  britt  le  fracture  models  even  though  such  models  have  been 
used  in  computations  for  many  years.  It  is  necessary  to  examine  this 
situation  further  to  resolve  these  inconsistencies.  It  may  be  possible 
to  extend  current  uniqueness-stability  theory  to  permit  some  form  of 
stress  discontinuity.  Another  possibility  is  that  acceptable  results  may 
be  obtainable,  for  some  types  of  ground  shock  problems,  if  a plasticity 
model  (such  as  contained  in  the  cap  model  routine  CAP75,  Ref.  [2]),  is 
used  to  model  tensile  failure.  These  ideas  and  others  can  be  examined  by 
means  of  a series  of  one-d imi'nsiona  1 wave  code  calculations  coupled  with 
Some  theoretical  studies  designed  to  focus  on  the  question  of  which  types 
of  tensile  models  are  appropriate  for  ground  shock  applications. 
Conclusions  and  recommc'ndat ions , based  on  the  above,  are  summarized  in 
Section  V. 
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II  BRITTLE  FRACTURE  MODEL  OF  SRI 


! 


i i 


« 


A.  Introduct  ton 

Although  the  brittle  fracture  model  of  SRI  was  originally  designed 
for  metals.  Refs.  [3,  4,  5],  it  has  been,  more  recently,  adapted  for 
ground  shock  applications  in  rocks.  Refs.  [6,  7,  8].  This  model  is  bused 
on  rate-dependent  concepts  of  fracture,  including  crack  nucleation  rates 
and  crack  growth  rates.  Stress-strain  behavior  of  tlie  material  is 
postulated  for  various  stages  of  damage,  including  full  fragmentation, 
defined  as  coalescence  of  cracks  .and  full  se|)aracion. 

A copy  of  the  model  routine,  BFRACT,  Ref.  [5],  w.is  obtained  and  has 
been  examined  theoretically,  as  well  as  computat  i on.a  1 1 y , for  ono- 
dimensional  problems.  The  b.asic  postulates  of  the  model  .and  the  results 
of  numerical  investigations  arc  discussed  in  this  section.  Some 
modifications  of  BFRACT,  based  (;n  the  assumption  of  sp.ati.al  isotropy,  are 
suggested . 

B.  Constitutive  ^(^i^atiqns  fjir  Cracked  _Materi_aJ 

The  BFRACT  model  routine  contains  many  elements  common  to  other 
constitutive  equation  routines  and  these  assumptions  are  summarized 
below . 

It  is  presumed  that  compressive  states  .are  represented  by  an 
appropriate  set  of  stress-strain  rel.ations  and  th.at,  at  some  point  based 
on  experimental  data,  the  tensile  behavior  recjuires  a representation  of 
cracks  and  their  propagation.  The  model  is  based  on  continuum  concepts 
and  tlius  it  is  assumed  th.at  the  microscopic  cracks  can  be  represented  by 
means  of  a few  simple  macroscopic  parameters. 
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liflien  fracture  begins,  the  specific  volume,  V,  at  i point  in  the 
continuum  is  given  by  the  sum 

V = V + V (]) 

s c 

where  is  the  specific  volume  of  the  solid  material,  and  is  tlie 

specific  volume  of  the  cracks.  The  pressure,  P,  at  a point  is  assumi-d 

to  be  related  to  the  pressure,  P , of  the  solid  material  by  the 

s 

Carroll-Holt  equation.  Ref.  1.9],  for  porous  materials,  i.e.. 


P V = P V 
s s 


(2) 


Since  tensile  strains  associated  with  brittle  fracture  are  initially 
small,  a linearized  Mie-Clrunelsen  equation  of  state  is  assumed  for  the 
solid  behavior 


P 

s 


•<0  V"  - V 

s s 


where  K is  the  ambient  bulk  modulus,  V the  initial  solid  specific 
o o 

volume,  r the  Gruneisen  ratio  and  E the  internal  energy.  Internal 
energy  changes  are  assumed  independent  of  the  "energy  of  distortion" 
so  that 


(3) 


dE  = - P dV 


(4) 


The  deviatoric  stress  tensor,  ^ point  is  calculated  by 

means  of 


ds 


ti 


= 2 u 


red 


(dc . . - 
ij 


dV 


) + 0, 


t.) 


(5) 


where  e . . 

1.1 


is  the  strain  tensor. 


S.  . 

i.I 


the  Kronecker  delta. 


6 . . is  a 
ij 
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r- 


correction  due  to  rotation  (large  displacement  theory)  and  ii  is  the 

red 


reduced  shear  modulus 


*) 


U = M (1  - 1.88  ■) 

red  o V 


*) 


(6) 


The  shear  strength  , Y,  is  assumc'd  to  decrease  linearly  with  increasing 
porosity  or  crack  volunii'  sc  tliat 

V 

Y = Y (1  - 4 ') 

o Y 

whore  a von  Mises  yield  loiuiilion 


\f7  = V.\/3 


is  assumed  and  .1^  is  tlu'  sti.ond  invariant  of  the  stress  devi;itors. 


It  remains  only  to  specify  an  equation  for  the  specific  volume  of 
the  cracks,  V . Vhir  a general  distribution  of  cracks  one  may  write 


V (k)  dN  (R) 
c c 


where  v is  the  volume  of  cracks  with  radius  R,  and  dN  is  the  number 
c c 

density  of  cracks  witli  radius  R.  The  Sneddon  re  1 at  ionsli ip  for  tlie 
elastic  opening  of  a penny-shaped  crack.  Ref.  [10],  is  assumed  so  tliat 
the  ellipsoidal  volume  of  the  crack  is 


_ 16  (1  - V^) 


(7) 


(8) 


(9) 


(10) 


*)  * 

The  factors  1.88  and  4 in  (J  , and  Y are,  in  fact,  material 

red 

parameters,  but  were  chosen  as  constants  in  all  SRI  work  so  far. 
See  Ref. [4,  p.  20]  for  discussion. 


7 


where  V is  Poisson's  ratio,  E is  Young's  modulus,  and  0 is  the  tensile 
stress  normal  to  the  plane  of  the  penny-shaped  crack. 

C.  Nucleation_  and  Growth  Rates 

The  laws  wiiich  govern  the  size  distribution  and  growtii  rates  of 
crack  radii  arc  the  relations  which  distinguish  the  SRI  model  from  other 
theories.  A distribution  function  for  the  Inherent  flaws  (penny-shaped 
cracks)  in  a brittle  substance,  as  well  as  the  laws  governing  their 
continued  nucleation  and  growth  under  tensile  stress,  are  postulated, 
witli  the  SRI  experimental  program  providing  the  material  parameters. 

The  usual  size  distribution  of  penny-shaped  cracks  seen  in  brittle 
materials  is  an  exponential  form 

= N(t)  exp  ( - R/R^)  (II) 

where  is  the  number  density  of  cracks  with  radii  greater  than  R,  N(t) 
is  the  total  number  density  at  time  t,  and  is  the  shape  parameter  of 
tiie  crack  size  distribution.  Tiie  crack  radius,  R,  is  also  a function  of 
time  which  will  subsequently  be  defined. 


The  rate  at  which  new  cracks  nucleate  under  the  influence  of  a 
*) 

tensile  stress  is  taken  to  be 


0-0 


N ex|> 
o 


0 < o 


N = 


(12) 


o > o 


'*y  ■ 

Stress  is  assumed  positive  in  compression. 
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where  N is  the  initial  nucleation  rate,  o is  the  nucleation  stress 
o no 

*) 

threshold  and  is  the  nucleation  stress  sensitivity  parameter. 

The  crack  growth  law  found  to  give  best  agreement  with  experimental 
data,  even  for  brittle  materials,  is  of  a viscous  form  (see  Ref.  [11]) 


T,  (a-a  )R  , o<a 

1 go  go 


R = 


(13) 


CJ  > o 


go 


where  T,  is  a constant  related  to  the  "viscosity"  and  o is  the  g.rowtli 
1 go 


threshold  stress. 


D.  FuJ  ^^Fragmejitaticm 


Up  to  this  point,  it  has  been  tacitly  assumed  tliat  the  growing 
cracks  do  not  influence  one  another.  To  go  from  mild  to  severe  damage, 
the  coalescence  of  cracks  and  the  resultant  formation  of  fragments  whiih 
cannot  support  tensile  stresses  must  be  accounted  for.  This,  is  done 
through  the  following  heuristic  argument. 


If  it  is  assumed  that  the  radius  of  a fragment,  R^.,  is 
proportional  to  the  size,  R,  of  the  cracks  forming  it,  and  tliat  the 


number  density  of  possible  fragments,  , is  proportional  to  the  number 

densltv  of  cracks,  N , two  material  constants  are  defined 

c 


(14) 


*) 

Tin-  apparent  d i st:ont  i nui  ty  in  the  nucleation  rate  function  is  not 

significant  since  a is  taken  to  be  the  tensile  stress  at  which 
no 

tile  crack  volume  first  becomes  non-zero. 
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and 


N 


(15) 


Experience  shows  y to  be  near  one  and  3,  reflecting  fragment  sliape,  to 
be  one-quarter  for  eight-sided  fragments,  as  are  most  often  assumed.  A 
third  parameter,  T^,  relating  relative  fragment  volume,  v^ , to  radius, 
is  also  of  order  unity. 


v^  = Tj,  R^' 


(16) 


Then,  the  relative  volume  of  all  fragment.?  is 


"f  = 


Vf  dN^ 


Upon  inserting  Eqs.  (14)  - (16), 


Vj.  = TpY  3 


dN 


(17) 


The  number  density  of  cracks  with  radius  R,  dN^ , is  obtained  from  Eq . 

(11). 

While  the  function  V^.  varies  continuously  frcjm  slightly  cracked  to 
fully  fragmented  material  states,  it  is  thought  to  overestimate 
fragmentation  in  the  early  stages  of  damage.  It  is  likely  that  some 
threshold  density  of  cracks  must  be  reached  before  the  likelihood  of 
cracks  intersecting  sufficiently  to  form  fragments  becomes  significant. 

This  threshold  value,  V^,  defining  the  beginning  of  crack  coa  1 e.scence , 
is  given  in  terms  of  the  material  constants  discussed  above  and  T^ , a 
measure  of  material  brittleness.  Refs.  [5,  11] 

T 

-/  (18) 
c 
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The  actual  relative  fragment  volume,  , is  defined  then  as 


Vf  - 

a ^ _f__  _* 

' 1 - v'. 


= 0 


V < V 
f - * 


(19) 


Total  porosity,  V , is  the  sum  of  the  crack  volumes  determined  from 

V 

Eqs.  (9)  - (13)  and  the  portion  of  the  solid  volume  which  has  fragmented 


V = V + V 
V c f s 


(20) 


(The  three  quantities,  V , V , V , are  specific  volumes.  V^‘  is  a 

V c s f 

re  1 a t i ve  vol ume ) . 


E.  BFRACT  Algor itlim 

While  the  nucleation  and  growth  laws  are  both  given  as  functitms  of 
the  normal  stress,  specifying  this  quantity  for  each  propagating  crack 
would  be  enormously  cumbersome  and  some  simplifying  assumptions  are  made. 
The  above  continuum  model  is  implemented  in  BERACT  by  first  assuming  tliat 
each  computat  ion.al  cell  can  be  broken  into  .an  array  of  crack  orientations 
or  "bins",  Ref.  [4,  p.  21].  For  two-dimensional  behavior,  each  bin  is 
defined  as  containing  cr.icks  perpendicular  to  (])'  (in  the  y plane), 
and  to  (perpendicular  to  the  x,  y plane),  where  the  index,  i,  refers 
to  the  i-th  bin.  Then,  direction-dependent  damage  functions  .are 
specified  relative  to  these  bin  orientations.  Thus,  a*",  is  the  tensile 
stress  in  the  direction  (i})^,  4>^)  t is  the  radius  of  a crack 
perpendicular  to  and  growing  under  the  influence  of  and  dN^  is  the 

time-dependent  number  density  of  cracks  with  radii  equal  to  R^. 


It  is  useful  to  see  the  transition  from  continuum  expression  to 
algorithm  for  a specific  function.  Integrating  Eq.  (13)  gives 


11 


rt 


R(t)  = R{o)  exp  [ 


T,  (a  (t')  - 0 ) dt' 

1 RO 


(21) 


for  the  radius  of  an  arbitrary  crack  under  the  influence  of  tensile 

stress  a (<  a ).  If  the  time  of  observation  is  taken  incrementally 
go 

close  to  a time,  t , at  which  the  radius  was  known  to  be  R , then 
o o 


R(t)  = R(t  + At)  - R exp  (T^  (a  - a ) At] 
o o 1 go 


(22) 


where  o is  the  average  stress  during  the  time  interval  At.  Using  the 
bin  breakup. 


R^(t)  ~ R ^ exp  [T  (o^,  - O ) At] 
o 1 (pw  go 


(23) 


Volvimetric  damage-related  functions  require  summation  over  the  bins 
in  each  cell.  So,  for  example,  the  relative  fragment  volume,  defined  in 
Eq . (17),  becomes 


3 ^ (R^)^  dN^, 


(24) 


F-  Storage  Requirements 


One  consequence  of  this  partitioning  into  bins  is  the  need  to  save 
at  least  "b"  memory  parameters  per  computational  element,  whore  b is  the 
number  of  bins.  In  fact,  two  time-dependent  functions,  N^(t)  and  R^(t), 
must  be  stored  for  each  bin.  In  addition,  the  void  volume,  V , and  the 

V 

unfragmented  volume  fraction,  f , are  stored,  where 

u 


f 

u 


1 


V 


a 

f 


making  a total  of  2b  + 2 memory  parameters  per  element.  As  BFRACT  is 
currently  implemented,  five  bins  (ler  element  are  assumed,  making  a total 
of  twelve  items  per  element  for  the  character i/.at  ion  of  brittle  tensile 


fracture. 
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G. 


Possible  Modifications  of  BFRACT 


The  function  of  BFRACT  is  to  calculate  the  time-dependent  growth  of 
damage  caused  by  the  application  of  tensile  stress,  while  accounting  for 
the  strength  reduction  caused  by  the  developing  damage.  The  routine 
functions  well  for  small-scale  shock  loading  problems  and  has  been  used 
successfully  to  model  plate  slap  experiments.  Ref.  [5],  armor 
penetration  problems.  Ref.  [5],  and  laboratory-dimension  crater 
detonations.  Ref.  [3]. 

However,  for  large-scale  applications,  such  as  many  ground  shock 
problems,  a scheme  which  requires  more  than  ten  memory  parameters  to 
model  tensile  behavior  appears  too  costly.  This  seems  especially  true 
when  the  uncertainties  in  modeling  the  in  situ  material  are  considered. 
For  this  reason,  it  was  decided  to  examine  ways  of  simplifying  BRFACT 
while  retaining  the  nucleation,  growth  and  fragmentation  concepts.  The 
hope  is  that  a significant  portion  of  the  modeling  capability  can  be 
preserved  but  at  a more  modest  cost  in  computer  storage. 


Two  approacnes  toward  simplification  could  be  considered,  depending 
on  the  type  of  fracture  occurring.  In  the  case  of  brittle  fracture 
(e.g.  of  rocks),  cracking  tends  to  be  highly  localized  and  iriented 
normal  to  the  axis  of  maximum  principal  stress.  In  contrast,  ^uc_H_l_e 
fracture  (e.g.  of  many  soils  or  metals)  occurs  in  a form  which  is 
characterized  by  spherical  voids  and  hence  is  associated  with  spatial 
isotropy . 


Thus , 

in  place  of 
criterion, 
s Ign i f leant 


for  brittle  fracture,  the  use  of  the  maximum  tensile 
would  seem  appropriate.  The  maximum  priP'  ioal 
if  implemented  witlmut  the  use  of  bins,  should  lead 
reduction  in  storage  reqtiirements . Similarly,  for 


St  ress 

st  ress 

to  a 
duct  i le 
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fracture,  the  use  of  a combination  of  stress  invariants  (e.g.  the 

I 

pressure  and  J^)  in  place  of  would  seem  appropriate  and,  if 

implemented  without  bins,  should  lead  to  a reduction  in  storage 
requirements. 

In  fact,  as  the  continuum  concepts  underlying  BFRACT  have  actually 
been  implemented  in  difference  form  bv  SRI,  a.,  is  variously 

4>4' 

approximated,  depending  on  its  sensitivity  in  a calculation  and  upon 
the  brittle  material  itself.  In  some  expressions  it  is  taken  as  a 
function  of  the  solid  pressure,  in  others  as  some  linear  combination  of 
the  solid  pressure  and  the  most  tensile  stress  deviator.  Thus, 
examining  different  approximations  for  o.,,  in  the  search  for  a more 
economical  BFRACT,  seems  qoite  appropriate. 


r 1 1 ON  E-  D I MENS  IONA  I.  STUDIES  BFRACT 

A.  Introduction 


In  the  preceding  section,  the  bnsic  continuum  relations  underlying  the 
model  routine  BFRACT  were  reviewed,  anil  some  possible  modifications  of  the 
algorithm  were  discussed.  In  this  section,  the  results  of  some  simple  one- 
dimensional  calculations  with  BFRACT  are  presented,  and  the  feasibility  of 
some  of  the  modifications  is  tested.  At  the  same  time,  the  rate-dependence 
of  the  model  routine,  and  its  sensitivity  to  parameter  changes,  is 
explored. 

B.  Armco_  I ron 

The  material  chosen  for  the  studies  is  Armco  iron,  n well  characterized 
brittle  metal  for  which  the  SRI  experimental  program  has  provided  a full 
set  of  model  parametirs.  The  reason  for  this  is  that  geological  materials, 
which  are  the  main  concern  i>f  the  present  study,  have  been  investigated 
less  fully  and  are  more  complex.  In  Table  1,  some  general  properties  of 
Armco  iron  .ire  listed.  In  Table  2.  those  par.imeters  necessary  for  the 
implementation  of  BFRACT  arc  given. 


Table  1.  Material  Properties  of  Armco  Iron  (From  Ref.  [5]) 


Mt'iss  density,  0 

o 

7.85  gm/ 

3 

cm 

Poisson's  ratio,  V 

0.2801 

Cruncisen  par.imeter,  V 

o 

1.690 

Bulk  modulus,  K 

o 

1.589  X 

in'^ 

dyn/cm^ 

Shear  modulus,  p 

o 

8.190  X 

dyn/ rm^ 

1 

Initial  vield  strength,  Y 

o 

5.  50  X 

10^ 

dyn/cm^ 

Bulk  sound  speed,  c 

o 

4.499  X 

10^ 

cm/ sec 

15 


Table  2. 


Fracture  and  Estimated  Fragmentation  Parameters  for  Armco  Iron 


Code 

Name  Tn 

Parameter 

Name 

Derivation 

Definition 

Value 

Units 

TSR(l) 

Growth  coefficient 

- 6.0  X 10~^ 

cm^/dyn/sec 

TSR(2) 

o 

go 

1 

(Growth  threshold 

- 2.0  X 10® 

dyn/cm^ 

TSR(3) 

Nucleation  crack 
size  distribution 

5.0  xio'^ 

cm 

TSR(4) 

N 

o 

Nucleation  rate 
coef  f icient 

12 

4.6  X 10 

3 

No. /cm  /sec 

TSR(3) 

o 

no 

Nucleation  threshold 

stress 

- 3.0  X 10^ 

dyn/cm^ 

TSR(6) 

Nucleation  stress 
sensitivity 

- 4.56  X 10^ 

dyn/ cm^ 

TSR(IO) 

6 

Ratio  of  the  number 

I of  fragments  to  the 
number  of  cracks 

0.25 

1 

! 

TSR(] 1) 

Y 

Ratio  of  fragment 
radius  to  crack 
radius 

1.0 

1 

TSR(12) 

V* 

Threshold  relative 
volume  of  potential 
fragments  to  produce 
actual  fragments 

0.2 

TSR(13) 

■^f 

Coefficient  relating 

1.0 



fragment  shape  to 
fragment  volume 

From  Kef.  [4,  p.  37 1 


Based  on  values  rliosen  for  armor  steel,  Ri’f.  (3,  p.  89] 


The  model  routine  BFRACT  was  driven  in  simple  uniaxial  extension  at 


a constant  strain-rate  of  1.8  x 10^  sec  Fig.  1 shows  a plot  of 


tensile  stress  versus  specific  volume  for  Armco  iron.  The  curve  labeled 
"SRI"  appeared  in  Ref.  [5,  p.  55],  and  the  curve  labeled  "WA"  was 
produced  at  Weidlinger  Associates  with  BFRACT,  using  the  parameters 
shown  in  Tables  1 and  2.  The  discrepancy  of  a few  percent  between  the 
two  curves  is  the  result  of  using  somewhat  different  parameter  values 
(due  in  part  to  rounded  values  given  in  reports).  Since  the  model  is 
quite  rate-sensitive  (see  bel<iw),  a few  percent  difference  is  not 
surprising. 


C. 1 Stress  Modifications  of  BFRACT 


As  a test  of  the  hypothesis  that,  f«>r  ground  shock  purposes,  BFRACT 
may  be  simplified  without  sacrificing  significant  modeling  capability, 
the  uniaxial  extension  test  was  rerun  with  the  following  modifications: 


Option  1 - Pressure  only 

Option  2 - Three  "bins"  instc.ad  of  five 


In  option  1,  the  quantity,  which  «'onlrols  the  development  of 

jjip 


tension-induced  f r.igmentat  ion , w.is  taken  to  be  independent  of  the  (most 
tensile)  deviatoric  stress.  Since  in  the  present  example  thi-  deviatoric 


1 


stresses  are  an  order  of  magnitude  smaller  than  the  pressure,  it  is  to 
be  expet'ted  that  setting  the  shear  modulus  to  zero  in  tensif’n  will  not 
affect  the  results  by  more  than  about  ten  per  cent. 

In  option  2,  the  number  of  bins  pc-r  elt-menl  was  reduced  from  five 
to  three  (the  minimum  number  which  can  be  accommodated  without  recoding 
or  suppressing  a stress  component).  It  is  expected  that  this  variation 
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FIG.  1 STRESS-V01.r!i;  PATH  FOR  STRAIN-RATE  UlADING  OF 

IRON  TO  fracme:ctatiox 


will  i>ri>duc,e  little  change  in  the  results  of  the  present  example. 

Fig.  2 confirms  these  expectations.  The  curve  labeled  SRI  is  the 
same  as  that  labeled  WA  in  Fig.  1.  The  three  other  curves  correspond  to 
option  1,  option  2,  and  options  J and  2 combined.  The  distinguishable 
quantity,  total  strain  to  full  fragmentation,  is  identical  to  within  four 
significant  figures,  while  peak  stresses  vary  by  less  than  ten  per  cent. 

C.2  Rate-Dependence  of  BFRACT 

BFRACT  is  a rate-dependent  model  devised  for  shock  loading 
applications  and  the  relationship  betv’een  strain-rate  and  damage  is 
highly  nonlinear.  For  illustrative  pvirposes,  the  strain-rate  in  the 
uniaxial  extension  problem  under  discussion  was  varied  over  several  orders 
of  magnitude.  Total  strains  to  full  fragmentation  and  peak  tensile 
stresses  are  listed  in  Table  3 and  their  variations  confirm  that  BFRACT 
is  strongly  rate-dependent. 

Table  3.  Rate-Dependence  of  BI'RACT  (Uniaxial  Extension) 


.Strain-R.ite 

Total  Strain 

Peak  Stress 

-1 

* 

(sec  ) 

(%) 

(kbars) 

5 

1.8  X 10 

6 

- 48 

1.8  X 10 ' 

2 

- 19 

1.8x10^ 

< 1 

- 8 

1 kbar  = . 1 GPa 


As  is  well  known,  many  metals  exhibit  ductile,  behavior  (prior  to 

fracturing  in  a brittle  mode)  under  quasi-static  testing  conditions,  as 

illustrated  in  Fig.  3,  Ref.  [3].  BFRACT  is  designed  to  model  brittle 

behavior  at  high  strain  rates,  but  not  the  ductile  transition  expected  in 

-3  -1 

the  quasi-static  regime  (E.  < 10  sec  ).  The  stress-strain  curves  of  tlie 
three  rapid  extension  tests  are  superimposed  on  Fig.  3 to  facilitate 
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FIG.  3 RATE  DEPENDENCE  OF  BFRACT  FOR  UNIAXIAI.  EXTENSION  OF  ARMCO  IRON 
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' comparison.  Again,  the  strong  rate-sensitivity  of  BFRACT  is  illustrated 

I by  the  dramatic  increase  in  peak  tensile  stress  with  strain-rate. 

One  of  the  potential  dangers  of  adopting  a model  such  as  BFRACT 

without  careful  examination  is  that  parameters  may  have  been  determined 
i 

i from  a set  of  experiments  with  a limited  range  of  strain-rates.  If  the 

I 

1 

! routine  is  used  indiscriminately  for  applications  (and  strain-rates) 

' other  than  those  for  which  it  was  envisioned,  difficulties  may  arise. 

I In  fact,  the  use  of  a rate-dependent  rodel  adds  another  dimension  to  the 

1 I modeling  problem  and  should  only  be  used  if  it  is  clearly  necessary.  An 

attempt  to  clarify  this  point  is  made  in  section  IV. 

The  one-dimensional  Lagrangi;in  wave  code  WONDY  (developed  at 
Sandia)  was  selected  as  a convenient  program  by  which  BFRACT  could  bt- 
examined  in  dynamic  calculations.  The  coding  and  data  structure  of 
WONDY  allow  for  the  inclusion  of  special  equation  of  state  and  material 
model  routines  with  a minimum  of  effort.  The  incompatibilities  which 
arise  between  BFRACT  and  WONDY  are  easily  overcome  for  the  present 
purposes. 

For  comparative  purposes,  tensile  behavior  was  modeled  in  three 

ways ; 


Case  1 ) BFRACT 

Case  2)  BFRACT  with  modification  option  1 (MODI) 

Case  3)  Instantaneous  tensile  fragmentation,  i.e.  the  stress 

instantaneously  reset  to  zero  when  the  tensile  strength 
of  a material  is  exceeded 
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Case  3)  represents  one  of  the  standard  tensile  fracture  options 
available  in  WONDY,  Ref.  [12].  When  the  stress  at  a cell  boundary  is 
reset  to  zero,  at  the  moment  of  fracture,  the  two  free  surfaces  created 
are  accommodated  by  the  introduction  of  an  additional  node. 

E . SRI  PI at e ^1  ap  Calculations 

In  the  1973  SRI  report.  Ref.  [5],  results  of  a tapered  flyer  plate 
experiment  - copper  impacting  steel  armor  - were  given  and  compared  to 
two-dimensional  code  calculations  in  order  to  test  the  predictive 
capability  of  BFRACT.  The  experimental  configuration  is  reproduced,  from 
Ref.  [5,  p.  60],  as  Fig.  4. 

For  the  coarse  computational  grid  initially  selected,  the  stress 
history  near  the  end  of  the  target  was  severely  eroded  and  a clear  spall 
signal  was  not  evident.  Therefore,  the  calculation  was  repeated  with  a 
finer  grid  in  a one-dimensional  calculation  in  which  only  the  central 
region  of  the  target  and  flyer  were  simulated.  The  discrepancy  between 
measured  and  computed  stress  histories  at  the  free  surface  of  the  target 
is  shown  in  Fig.  5,  Ref.  [5,  p.  64].  It  was  felt  that  agreement  between 
the  two  waveforms  might  be  improved  by  altering  the  material  parameters 
in  such  a way  as  to  del;iy  the  damage  process.  In  particular,  it  was 
suggested.  Ref.  [5,  p.  65],  that  the  spall  arrival  time  might  be 
critically  dependent  on  V^,  the  threshold  volume  (relative)  for 
fragmentation.  This  hypotliesis  was  not  tested  for  armor  steel,  but  such 
parameter  sensitivity  was  examined  for  Armco  iron,  as  described  in  the 
next  section. 

F . M a_^e  Slap  Pa r^wie ter  Sensi  t ijy  ity^  S tj.idjr 

In  the  performance  of  the  parameter  sensitivity  study,  it  would  have 


22 


been  useful,  for  comparative  purposes,  to  have  utilized  the  plate  slap 
configuration  of  Fig.  4.  However,  the  BFRACT  parameters  appropriate  to 


copper  were  not  directly  available  (though  they  could  have  been  obtained 
from  SRI)  and  copper  and  steel  are  less  rock-like  in  their  responses  to 
tension  than  the  highly  brittle  Armco  iron  previously  selected  for 
examination.  Thus,  one-dimensional  wave  code  calculations  were  made,  as 
part  of  the  present  study,  for  the  case  of  an  Armco  flyer  (of  0.1270  cm) 
impacting  on  an  Armco  target  (of  0.4233  cm)  with  a velocity  of 

4 

5.0  X 10  cm/sec.  The  size  and  time  scales  are  similar  to  the  SRI  one- 
dimensional  simulation  discussed  in  the  previous  section. 

The  free  surface  velocity  histories  of  the  Armco  target  are  shown 
in  Fig.  6 for  tlie  three  cases  of  tensile  modeling  listed  in  Section  D. 
For  this  simple  problem,  the  first  two  cases  - BFRACT  and  BFRACT  MODI  - 
are  practically  indistinguishable.  As  expected,  instantaneous 
fragmentation  results  in  a slightly  earlier  s)iall  signal  arrival  time. 
The  velocity  histories  subsequent  to  the  spall  time  are  similar  but  the 
instantaneous  case  results  in  a slightly  higher  average  velocity.  This 
is  probably  due  to  the  fact  that  free  surface  velocity  continues  to 
decrease,  once  the  loading  wave  reflect.s,  until  the  arrival  of  the  spall 
signal.  Thus,  the  earlier  t!\o  time  <if  spall,  the  less  the  velocity  has 
decreased . 

The  sensitivity  of  BFRACT  to  rate-dcpendc'iice  was  discussed  earlier 
and  is  quite  strong,  as  shown  in  Table  3 and  Fig.  3.  In  order  to  study 
material  parameter  sensitivity,  the  parameter  list  can  be  separated  into 
a fracture  set: 

[TSR(J),  .1  = I,  61 
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6 COMPARISON  OF  COMPUT'^D  FRCi;  SIRFACF.  VFFOCITY  HISTORIES  OF  ARMCO 
TARCET  FOR  THREE  CASES  OF  TENSILE  MODELINC; 


Jo.i 

'(12)=<0.2 


COMI’ARISON  OF  COMPETED  FREE  SURFACE  VELOCITY  HISIORIES  OF  ARMCO 
TARCET  FOR  DIFFERINC  VALUES  OF  PARAMETER  1(12) 


and  a fragmentation  set: 


[TSR(J),  J = 10,  13] 

(See  Table  2 for  the  physical  significance  of  these  parameters.) 

The  value  of  V^,  [T(12)],  was  varied  from  0.1  to  0.5,  the  original 
choice  being  0.2.  The  free  surface  velocity  histories  are  shown  in 
Fig.  7 for  these  three  values  of  V^.  No  change  was  detected  in  spall 
signal  arrival  time  and  virtually  no  change  occurred  in  the  post-spall 
signal,  contrary  to  expectations. 

In  an  effort  to  detect  sensitivity  to  the  fragmentation  parameter 
T^,  [T(13)]  was  Increased  from  1.0  to  4.0  (the  latter  being  an 
alternate  choice  in  various  SRI  BFRACT  calculations.)  The  post-spall 
velocity  history,  shown  in  Fig.  8,  became  slightly  oscillatory  but  the 
difference  is  not  considered  significant  in  the  context  of  ground  shock. 

As  another  test  of  sensitivity,  some  of  the  fracture  parameters 
were  varied  dramatically,  being  assigned  values  appropriate  to  armor 
steel  rather  than  to  Armco  iron.  As  shown  in  Table  4,  parameter 
values  changed  by  as  much  as  four  orders  of  magnitude.  While  the  spall 
signal  arrival  time,  from  Fig.  9,  became  less  clearly  defined,  the  post- 
spall signal  was  not  altered  very  much. 

Thus,  the  picture  that  emerges  from  this  limited  study  is  of  the 
BFRACT  model  being  highly  rate-sensitive,  but  much  less  sensitive  to  a 
number  of  the  material  parameters. 
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FK;.  8 COMPARISON  OF  COMI’UTED  FREE  SURFACE  VELOCITY  HISTORIES  OF  ARMCO 
TARGET  FOR  DIFFERING  VALUES  OF  PARAMETER  T(]3) 


FIG.  9 


COMPARISON  OF  COin’UTEI)  FREE  SURFACE  VELOCITY  HISTORIES  OF  ARMCO 
TARGET  FOR  TWO  SETS  OF  FRACTURE  PARj\METERS  T(I)  - T(6) 
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Table  4. 


Dynamic  Fracture  Parameters  for  Armor  Steel 


*) 


and  Arraco  Iron 


Code 

Dos i^nat ion 

Name  In 

Derivations 

Units 

Armor 

Steel 

Armi-o 

I ron 

TSR(l) 

Tj  , Crowth 

coef  f icient 

cm^/dyn/sec 

- 0.00012 

- 0.0006 

TSR(2) 

0 or  P 

Pd'  go 

Growth  threshold 

dyn/ cm^ 

0.0 

, -3 

- 2.0  X 10® 

TSR(3) 

R^,  Nucleation 
size  parameter 

cm 

4.0  X 10 

5.0  X 10  ^ 

TSR(4) 

N , Threshold 
o 

nucleation  rate 

No. cm  /sec 

5.0  X 10® 

12 

4.6  X 10 

TSR(5) 

0 or  P , 

no  no 

Nucleation 

threshold 

2 

dyn/cm^ 

- 1.12  X 10^° 

- 3.0  X 10*^ 

TSR(6) 

Nucleation 
sens  it ivi ty 

dyn/cm 

- 7.4  X 10® 

- 4. 56  X 10^ 

*) 

From  Ref.  [5,  p.  89] 
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experimental  tests,  has  been  examined  in  some  detail  as  a candidate  model 
for  tensile  behavior  of  geological  materials  in  ground  shock  calculations. 
It  is  important  at  this  point  to  try  to  specify,  us  clearly  as  |)ossible, 
the  type  and  scale  of  problems  which  are  intended  by  the  phrase  "grt)und 
shock  applications."  This  is  done  in  the  present  section  and  it  is  argued 
that,  for  practical  reasons,  tensile  behavior  should  be  modeled  as  a rate- 
independent  effect  for  this  class  of  problems. 

B.  l^aijj^ejScale  Grj)und_^h_ock  Pjoblems 

There  are  many  land-based,  near  surface,  defensive  systems  for  whiih 
ground  motion  effects  from  large-scale  explosions  are  of  concern  (see 
Fig.  10).  Such  systems  would  not  be  expected  to  survive  a one  megaton 
nuclear  surface  burst  at  a distance  of  less  than,  say,  200  meters  because 
this  distance  is  less  than  tlie  crater  radius  in  ordinary  geologies. 
However,  beyond  two  or  three  crater  radii,  it  is  possible  to  design 
protective  structures,  if  the  exjiected  airblast  and  ground  shock  effects 
can  be  assessed.  It  is  a currently  accepted  principle  that  the  cost  of 
such  systems  are  kept  to  a more  reasonable  level  if  site-specific,  tlireat- 
specific  case  studies  are  made.  In  order  to  reduce  the  number  of  case 
studies  to  a tolerable  level,  an  atti'mjit  must  be  made  to  characterize  the 
most  critical  parameters  for  such  problems. 

Tbe  threat-specific  parameters  are,  I rom  tlie  standpoint  of  round 
shock,  relatively  few  in  number  compared  to  tbe  site-specific  parameters. 
For  many  ap])  I icat i ons , the  size  and  location  of  the  explosion  are  all 
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TYPICAL  RANGES  OF  INTEREST  (1  MT) : 0.3  to  3 kilometers 

TYPIGAL  TIMES  OF  INTEREST  (I  MT)  : 0.1  to  10  seconds 


FIG.  10  LARGE-SCALE  GROUND  SHOCK  PROBLEMS 
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that  are  required.  (Otlier.  roustily  equivalent,  parameters  are  also  used.) 
Given  these  [larameters,  it  is  presumc‘d  tliat  an  air  overpressure  time 
history  whiel)  is  independent  of  azimuthal  direction  can  be  defined  and 
used  as  an  applied  surface  loadinj;.  This  procedure  has  many  hidden 
assumptions  such  as: 

a)  details  of  the  explosive  process  are  not  important, 

b)  details  of  crater  formation  are  not  Important, 

c)  details  of  air-ground  interactions  are  not  important. 

The  site-specific  parameters  are  usually  characterized  by  a site 
profile,  which  is  a list  of  geological  materials  versus  depth  from  ground 
surface.  Implicit  in  this  procedure  are  hidden  assumptions  such  as: 

a)  the  site  can  be  considered  to  be  horizontally  layered, 

b)  each  layer  is  homogeneous  and  isotropic  (perhaps, 
transversely)  , 

c)  a material  model  for  each  layer  can  be  determined. 

Typically,  with  current  continuum  mechanics  assumptions,  the  number  of 
materials  is  between  two  and  twelve,  and  the  number  of  model  parameters 
for  each  material  is  three  or  mc're.  Thus,  the  number  of  model  parameters 
needed  to  define  a typical  ground  shock  calculation  is  between  six  and  .1 
few  hundred.  If  the  assumjit  ion  of  horizontal  layering  is  considered 
invalid  or  the  models  are  made  quite  complicated,  then  the  number  of 
parameters  can  increase  even  furtlier.  Such  complicated  calculations  are 
sometimes  made  for  research  purposes  but  should  be  avoided  unless  shown 
to  be  absolutely  necessary  since  the  effort,  both  experimentally  and 
computationally,  also  will  Increase  considerably.  Since  the  introduction 
of  ratc-dopendonce  will  add  a new  dimension  to  the  modeling  process,  it 
is  worthwhile  avoiding  it,  if  possible,  in  ground  shock  applications. 
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C.  Time  Scales 

The  time  scale  associated  with  the  rate-dependent  tensile  behavior 
of  Armco  iron  can  be  estimated  by  studying  the  results  of  Sections  II 
and  III.  Fig.  4 shows  that  the  spall  signal  arrival  times  for  the 
instantaneous  and  rate-dependent  models  differ  by  approximately  0.1 
microseconds.  The  behavior  of  the  iron  in  this  problem  can  be  described 

4 

as  follows.  The  flyer  plate,  Impacting  at  a speed  of  5 x 10  cm/sec, 
induces  a compressive  stress  pulse  with  an  amplitude  of  more  than  90 
kbars  (9  GPa)  in  the  iron.  The  stress  pulse  travels  to  the  free  end  of 
the  target,  attenuating  to  about  60  kbars  (6  GPa),  reflects  as  a 
tensile  wave,  and  creates  tensile  stresses  (locally)  of  about  50  kbars 
(5  GPa).  These  tensile  stresses  can  be  sustained  for  only  a short  time 
however.  As  shown  in  Table  2 (Section  III),  for  uniaxial  extension,  a 
peak  tensile  stress  of  48  kbars  (4.8  GPa)  corresponds  to  a strain-rate 
e = 1.8  X 10^  sec  \ and  a total  fragmentation  strain  Cj-  = 0.06.  At 
this  strain-rate,  the  target  will  fragment  in  a time 

t^  = — = 0.33  X 10  sec. 

F. 

Thus,  the  stress  can  be  expected  to  reach  its  most  tensile  value  of  50 
kbars  (5  GPa)  in  about  half  this  time  .ind  then  relax  back  to  zero  (see 
Fig.  1).  It  is  quite  reasonable  then  that  the  strain-rate  effect  leads 
to  a difference  of  0.1  microseconds  in  spall  signal  arrival  times. 

Using  this  case  as  a guide,  it  can  be  seen  that  a flyer  plate  of 
lower  velocity  will  induce  a lower  strain-rate  and  the  time  delay  will 
increase  somewhat.  (This  trend  has  been  cibserved  bv  performing 
calculations  with  flyer  plate  speeds  of  2.5  x IP"*  cm/sec  and 

4 

1 X 10  cm/sec.)  However,  since  the  total  i r.igmentat ion  strain  also 
decreases,  the  delay  times  do  not  change  vis  rapidly  as  the  strain-rate. 
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For  example,  if  the  induced  tensile  stress  is  about  8 kbars  (.8  CPa) 

Table  2 shows  that  the  fragmentation  time,  as  estimated  above,  will  be 

about  5 microseconds,  and  therefore  a spall  signal  delay  of  ^ or  1 

microseconds,  as  compared  to  instantaneous  fracture,  would  be  expected. 

Delay  times  much  longer  than  a few  microseconds  are  not  detectable  since 

the  nucleation  threshold,  0 , is  about  3 kbars  (.3  GPa). 

no 

Up  to  this  point,  only  Armco  iron  was  examined  because,  as  noted 
earlier,  there  is  relatively  little  rate-dependent  data  available  on 
geological  materials.  Nonetheless,  it  is  possible  to  make  some  general 
statements  about  soils  and  rocks  for  use  in  ground  shock  problems.  As 
indicated  in  Ref.  [1],  soils  close  to  the  surface  have  tensile  strengths 
which  are  practically  zero  (less  than  1 bar)  while  competent  rocks  at 
depth  may  have  tensile  strengths  of  more  than  300  bars  (.03  GPa). 
Compressive  strengths  are  about  10  times  the  tensile  strength,  and 
modulus  values  are  100  to  1,000  tmes  the  compressive  strengths. 

Table  5 shows  representative  values  for  several  geological  materials. 
(See,  for  example.  Refs.  [13,  14]). 

Table  5 "Representative"  Properties  of  Geological  Materials 


Material 

Tensile 

Strength 

(kbar)* 

. . . 

Compressive 

Strength 

(kbar) 

Loading 

Modulus 

(kbar) 

"Soil" 

.0001  to  .001 

.001  to  .01 

.3  to  3 

Sandstone 

.01 

. 1 

30 

Shale 

.03 

. 3 

100 

Limestone 

. 1 

1 

300 

Slate 

. 1 

1 

300 

Granite 

. 2 

1 . 5 

500 

* 


kbar  = . 1 GPa 
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Tlu'st'  prDpt’rt  ies  Ciin  vary  considerably  and  are  presented  only  as  a rougli 
indii-ation  of  parameter  values  for  some  "typical"  materials. 

Data  on  rate-dependent  tensile  behavior  of  soils  and  rocks  is  qnitt- 
meager  but  what  little  there  is  indicates  that  these  materials  di>  luu 
require  rati--dependent  representation  in  large-scale  ground  shock 
calculations  of  the  type  depicted  in  Fig.  10.  This  is  indic.ited  by  the 
fact  that  time  scales  of  less  th.in  10  microseconds  are  needed  to  risolve 
rate-dependence  in  these  materials,  yet  ground  shock  calculations  are 
conducted  with  time  steps  one  or  two  orders  of  magnitude  l.irger.  TIu' 
time  stop  is  controlled  by  wave  prop.igation  effects  in  .1  grid  with  boxes 
of  the  order  of  one  meter  or  more  and  with  wave  speeds  of  six  kilometers 
per  second  or  less.  Time  steps  greater  than  100  microseconds  are  thus 
typical  in  large-scale  ground  shock  applications.  (On  the  other  hand, 
sm.ill-scale  problems,  e.g.  laboratory  specimens  of  a few  centimeters, 
require  microsecond  scale  time  steps  and  rate-d('pendent  effects  are 
detectable. ) 

I).  Stability  and  Unijgiieness 

In  view  of  the  above  time-scale  discussion,  it  would  be  preferable 
to  model  tensile  failure  by  means  of  a rate-independent  continuum 
description.  However,  it  is  important  that  certain  theoretic-al  conditions 
be  satisfied  to  guarantee  that  the  continuum  model  yield  a unique 
solution.  To  do  so  may  require  some  modeling  approximations  which  then 
need  examination  to  determine  their  range  of  validity. 

The  key  point  to  observe  is  that  tensile  failure  is  always  associated 
with  a phenomenon  which  can  be  characterized  as  brittle.  There  may  be 
various  ductile-like  transitions  but  ultimately  there  appears  a rapid 
relaxation  from  a measurable  non-zero  stress  level  to  a zero  stress.  This 
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relaxation  occurs  in  a time  frame  which  is  practically  instantaneous.  At 
this  point,  separation  of  the-  material  occurs  and  it  is  no  longer  a true 
"continuum".  However,  in  light  of  the  variability  of  geological 
materials,  it  is  unreasonable  to  try  to  keep  track  of  the  various 
particles  which  are  created  by  this  process. 

On  the  other  hand,  the  stability-uniqueness  theory  for  rate- 
independent  continua  is  Incompatible  with  the  instantaneous  stress  jumps 
(discontinuities)  needed  to  represent  brittle  tensile  failure  . One 
approach  to  this  situation  is  to  examine,  the  uniqueness  theory  and  try  to 
extend  the  theory  somewhat  to  permit  these  discontinuities.  This  is 
precisely  what  was  done  in  the  case  of  shock  waves  in  continua.  Ref.  [15, 
p.  490].  Shock  waves  are  stress  jumps  in  compression  (or  shear)  and  can 
be  understood  physically  as  the  limit  of  viscous  material  behavior  as  the 
viscosity  becomes  extremely  small.  Mathematically  tliese  discontinuous 
functions,  called  weak  solutions,  are  obtained  as  the  limit  of  continuous 
solutions.  Numerically  the  solutions  are  discrete  functions  whicti 
approximate  the  weak  solutions;  viscosity  is  sometimes  included  in  the 
scheme  in  order  to  maintain  numerical  stability.  For  linear  equations 
the  theory  is  ci)mplete,  and  for  nonlinear  problems  a number  of  special 
cases  (mainly  one-dimensional)  have  been  worked  out  completely  so  that  a 
phys ica 1 1 y,  mathemat ic;i  I 1 y and  numerically  coherent  theory  exists  for  a 
large  class  of  problems.  However,  brittle  tensile  failure  is  quite 
different  from  a shock  wave,  and  the  theory  referred  to  cannot  simply  be 
transferred  to  this  situation. 

*) 

Private  communication,  l.S.  Sandler,  Weidlinger  Associates,  1976. 


35 


Lt  is  interesting  tlien  that  numerical  solutions  have  been  obtained 


using  models  which  represent  idealized  brittle  failure.  For  example,  the 
TENSOR  code.  Ref.  [16,  p.  198],  has  allowed  a brittle  fracture  option  for 
many  years  and  (presumably)  no  instability  has  been  traced  to  the  use  of 
tills  model.  As  another  example.  Ref.  [17]  presents  results  for  one- 

dimensional  spherical  wave  propagation  in  a material  which  is  linearly 

* ) 

elastic,  except  that  an  Instantaneous  brittle  fracture  model  is 
introduced  when  a specified  tensile  hoop  stress  develops.  There  is  no 
theoretical  basis  for  accepting  the  results,  but  since  they  appear 
physically  sensible  and  numerically  stable  there  may  be  a mathematical 
basis  for  using  certain  types  of  brittle  fracture  models.  One  possibility 
then  is  to  examine  this  numeric:al  model  to  sec  whether  the  solution 
corresponds  to  a mathematically  well-posed  jiroblcm  which  is  physically 
acce]itab  1 e . 

Another  approach  to  the  uniqueness  problem  is  to  try  to  model  brittle 
failure  by  means  of  plasticity  for  which  there  is  an  established 
uniqueness  tlieory.  The  liope  here  is  that  this  approx  i m.it  i on  (i.e.  ductili' 
failure  instead  of  brittle)  can  be  shown  to  be  accept. ible  for  usi-  in 
large-scale  ground  shock  applications.  This  is  thi-  ti'chnitpit'  which  is 
used  in  the  current  cap  model  routine.  Ref.  [2].  A siTies  of  simple  I'lie- 
dimensional  calculations  can  be  made  wliicli  will  allow  this  mi>deling 
technique  to  be  examined.  The  stability  and  uni(|u(.'ness  thi’orv  for  hrittli' 
fracture  can  be  studied  at  the  s.ami'  time  (piite  convenientlv  and 
recommendations  concerning  these  ideas  .arc  discussed  in  Scctic'ii  V. 


*) 


The  brittle  fracture  model  in  T'ENSdR  reduces  to  the  model 
Ref.  [15]  for  the  case  of  spherical  symmetry. 


used  in 
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V CONCJr-PALON 

A.  Sunima  ry 

As  recommended  in  Ref.  [1],  tlie  SRI  rate-dependent  brittle  fracture 
model  (as  implemented  in  BFRACT)  has  been  examined  by  means  of  some  simple 
one-dimensional  calculations  using  parameters  for  Armco  iron.  The  model 
is  strongly  rate-dependent  but  time  scales  of  a few  microseconds  are 
needed  to  resolve  these  rate  effects.  From  the  meager  data  available  on 
geological  materials,  rate  effects  in  rocks  are  believed  to  be  of  the  same 
time  scale.  Since  large-scale  ground  shock  calculations  are  typically 
performed  with  time  steps  larger  than  100  microseconds,  it  is  impractical 
to  consider  using  such  a model  for  these  problems  and  other  modeling 
techniques  should  be  considered. 

Ultimately  tensile  failure  always  apiiears  as  a brittle  phenomenon 
wherein  a practically  instantaneous  stress  drop  occurs.  In  its  idealized 
rate- i ndependent  form,  brittle  fracture  is  incompatible  with  stability  and 
uniqueness  theory  for  a continuum.  (Once  fracture  occurs  tlie  material  is 
no  longer  even  a continuum.)  However,  two  [)ossib!e  approaclies  sliould  be 
examined  since  they  may  be  acceptable  for  brittle  tensile  modeling  in 
large-scale  ground  shock  applications.  Fortunately,  these  two  ap|ironches 
can  be  studied  at  the  same  time  quite  conveniently. 

B.  f^(^ommendii_tJ_on_s 

It  is  recommended  that  a series  of  one-dimensional  problems  be 
examined  in  order 


(i)  to  determine  if  it  is  possible  to  use  an  ideal  izi'd  brittle 

fracture  model  in  a wave  code  with  the  guarantee  that  a well- 
posed  problem  (with  a unif|ue  solution)  is  being  solved. 


(ii)  to  test  tile  raiig,e  of  validity  of  using  a plasticity  model  in 
place  of  a brittle  fracture  model  in  ground  shock 
app] i cat  ions. 

One-dimensional  wave  propagation  will  be  studied  with  a rate- independent 
brittle  fracture  model.  This  numerical  model  will  be  examined 
theoretically  to  try  to  determine  if  it  re|>resents  a iihysically 
reasonable,  mathematically  well-posed  problem.  The  plasticity  model 
which  is  currently  used  in  the  cap  model  routine.  Ref.  [2],  will  also  be 
used  in  the  same  wave  code  and  the  results  will  be  examined  with  large- 
scale  ground  shock  applications  in  mind.  By  this  means  it  is  hoped  that 
an  acceptable  rate-independent  tensile  model  can  be  developed. 
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AF  Institute  of  Technology,  AU 

ATTN;  Library  AFIT,  Bldg.  640,  Area  B 

AF  Weapons  Laboratory,  AFSC 

ATTN:  DES-S,  M.  A.  Plamondon 
ATTN:  SUL 

ATTN:  DES,  G.  Wayne  Ullrich 
Headquarters 

Air  Force  Systems  Command 
ATTN:  DLCAW 

Commander 

Foreign  Technology  Division,  AFSC 


ATTN: 

NICD  Library 

Hq.  USAF/IN 

ATTN: 

INATA 

Hq.  USAF/PR 

ATTN: 

PRE 

Hq.  i:SAF/RD 

ATTN; 

RDQSM 

Commander  Rome  Air  Development  Center,  AFSC 
ATTN:  EMTLD,  Doc.  Lib. 

SAMSO/MN 

ATTN : MNN 

Commander  in  Chief 
Strategic  Air  Command 

ATTN:  NRT,  STINFO  Lib. 

DEPARTMENT  OF  ENERGY 

Department  of  Energy 
Albuquerque  Operations  Office 

ATTN:  Doc.  Con.  for  Tech.  Library 

Department  of  Energy 

Division  of  Headquarters  Services 

ATTN:  Doc.  Con.  for  Class.  Tech.  Lib. 


Commander 

Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 

ATTN:  Tech.  Lib. 

Commanding  Officer 

Naval  Underwater  Systems  Center 

ATTN:  Code  EM,  Jack  Kalinowskl 


Department  of  Energy 
Nevada  Operations  Office 

ATTN:  Doc.  Con.  for  Tech.  Lib. 

University  of  California 
Uiwrence  Livermore  Laboratory 

ATTN:  Larry  W.  Woodruff,  L-96 
ATTN:  Tech.  Info.,  Dept.  L-3 


President 

Naval  War  College 

ATTN:  Tech.  Library 


Los  Alamos  S'^ientlfic  Laboratory 

ATTN:  Doc,  Con.  for  Reports  Lib. 
ATTN:  Doc.  Con.  for  R.  J.  Bridwell 


Commanding  Officer 
Naval  Weapons  Evaluation  Facility 
ATTN:  Tech.  Lib. 


Sand i a Laborator ics 
Livermore  Labontory 

ATTN:  Doc,  Con.  foi  Tech.  Lihrarv 


Director 

Strategic  Systems  Project  Office 
ATTN:  NSP-43,  Tech.  Lib. 


Sandia  Laboratories 

ATiN:  Doc.  Con.  for  3141,  Sandia  Rpt . Col). 

ATTN:  L.  Hill 
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Union  Carbide  Corporation 
Holifieid  National  Laboratory 

ATTN:  Civil  Def.  Res.  Proj. 
ATTN:  Doc.  Con,  for  Tech.  Lib. 

OTHER  GOVERNMENT  AGENCY 

Department  of  the  Interior 
Bureau  of  Mines 

ATTN:  Tech.  Lib. 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corporation 

ATTN:  Tech.  Info.  Services 


Institute  for  Defense  Analyses 

ATTN:  IDA  Librarian,  Ruth  S.  Smith 

Kaman  AviDyne 

Division  of  Kaman  Sciences  Corp. 

ATTN:  Tech,  Lib. 

ATTN:  E.  S.  Criscione 
ATTN:  Norman  P.  Hobbs  j 

Kaman  Sciences  Corporation 
ATTN:  Library 

Lockheed  Missiles  & Space  Co.,  Inc. 

ATTN:  Tech.  Info.  Ctr.,  D/Coll. 

ATTN:  Tom  Geers,  D/52-33,  Bldg.  205 


Agbabian  Associates 

ATTN:  M.  Agbabian 

Applied  Theory,  Inc. 

2 cy  ATTN:  John  C.  Trulio 

Avco  Research  & Systems  Group 

ATTN:  Research  Lib.,  A830,  Rm.  7201 


Lovelace  Foundation  for  Medical 
Education  and  Research 

ATTN:  Asst.  Dir.  of  Res.,  Robert  K.  Jones 
ATTN:  Tech.  Lib. 

McDonnell  Douglas  Corporation 
ATTN:  Robert  W.  Halprin 

Merritt  CASES,  Inc. 


Battelle  Memorial  Institute 

ATTN 

J . L . Merritt 

ATTN:  Tech.  Lib. 

ATTN 

Tech.  Lib. 

The  BDM  Corporation 

Nathan  M.  Newmark 

ATTN:  Tech.  Lib. 

Consulting 

Engineering  Services 

ATTN 

Nathan  M.  Newmark 

The  Boeing  Company 

ATTN:  Aerospace  Library 

Physics  International  Company 

ATTN 

Doc.  Con.  for  Tech.  Lib. 

California  Research  & Technology,  Inc. 

ATTN 

Doc.  Con.  for  E.  T.  Moore 

ATTN:  Ken  Kreyenhagen 

ATTN 

Doc.  Con.  for  Robert  Swift 

ATTN:  Tech.  Lib. 

ATTN 

Doc.  Con.  for  Larry  A.  Behrmann 

ATTN:  Sheldon  Shuster 

ATTN 

Doc.  Con.  for  Fred  M.  Sauer 

ATTN 

Doc.  Con.  for  Dennis  Orphal 

Calspan  Corporation 

ATTN:  Tech.  Lib. 

R&D  Associates 

ATTN 

Harold  L.  Erode 

Civil/Nuclear  Systems  Corp, 

ATTN 

Albert  L.  Latter 

ATTN:  Robert  Crawford 

ATTN 

Tech.  Lib. 

ATTN 

Robert  Port 

University  of  Dayton 

ATTN 

Henry  Cooper 

Industrial  Security  Supt-r.,  K[,-')n5 

ATTN 

Jerry  Carpenter 

ATTN:  Ha  I lock  F.  Swift 

ATTN 

William  B.  Wright.  Jr. 

ATTN 

J,  C.  Lewis 

University  of  Denver 

Co  1 ora  do  S em  i n.i  r y 

Science  Appl  icat  i(’ns,  Inc. 

Denver  Research  InstUuti- 

ATTN 

D.  E.  Maxwell 

ATTN:  Sec.  Officer  for  J.  Wisotski 

ATTN 

David  Bernstein 

EG&G,  Inc. 

Science  Applications,  Inc. 

Albuquerque  Division 

ATTN 

Tech.  Lib. 

ATTN:  Tech.  I.Ib. 

Southwest  Research  Inst  i tute 

Gard,  Inc. 

ATTN 

Wilfred  E.  Baker 

ATTN:  G.  L.  Neidhardt 

ATTN 

A.  R.  Won7.el 

General  Electric  Companv 

SRI  International 

TF.MPO-Center  for  Advanced  Studies 

ATTN 

George  R.  Abrahamson 

ATTN:  D AST AC 

ATTN 

Burt  R.  Gasten 

ATTN 

\ . Gvjpt.a 

ITT  Research  Institute 

ATTN:  R.  E.  Welch 
ATTN:  Milton  R.  Johnson 
ATTN:  Tech.  Mb. 
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Systems,  Science  & Software,  Inc. 

The  Eric  H 

Wang 

ATTN: 

Tech.  Lib. 

Civil  Engineering  Rsch,  Fac . 

ATTN: 

Donald  R.  Grine 

ATTN 

Neal  Baum 

ATTN: 

Ted  Cherry 

ATTN 

Larry  Bickle 

ATTN: 

Thomas  D.  Riney 

Weidl inger 

Assoc.  Consulting 

Engineers 

Terra  Tek, 

tnc . 

ATTN 

Melvin  L.  Baron 

ATTN: 

Tech.  Lib. 

ATTN 

Joe  Wright 

ATTN: 

Sidney  Green 

ATTN 

Lorraine  Whitman 

ATTN: 

Ahmed  Abou-Soyed 

Weigl inger 

Assoc.  Consulting 

Engineers 

Tetra  Tech. 

, Inc . 

ATTN 

J.  Isenberg 

ATTN: 

Li-San  Hwang 

ATTN: 

Tech.  Lib. 

Westinghouse  Electric  Corp. 

Marine  Division 

TRW  Defense 

& Space  Sys.  Group 

ATTN 

W . A.  \'o  1 7. 

ATTN: 

Tech.  Info.  Center,  5-1930 

ATTN: 

Pravin  Bhutta,  Rl-1104 

2 cy  ATTN: 

Peter  K.  Dai,  Rl-2170 

TRW  Defense 

& Space  Sys,  Group 

ATTN: 

E.  Y,  Wong,  527/712 

Universal  Analytics,  Inc, 
ATTN:  E.  I.  Field 
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